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Abstract 

We make a detailed analysis of both phenomenological and analytic back- 
ground for the "Brownian recoil principle" hypothesis (Phys. Rev. A 46, 
(1992), 4634). A corresponding theory of the isothermal Brownian motion 
of particle ensembles (Smoluchowski diffusion process approximation), gives 
account of the environmental recoil effects due to locally induced tiny heat 
flows. By means of local expectation values we elevate the individually neg- 
ligible phenomena to a non-negligible (accumulated) recoil effect on the en- 
semble average. The main technical input is a consequent exploitation of 
the Hamilton-Jacobi equation as a natural substitute for the local momen- 
tum conservation law. Together with the continuity equation (alternatively, 
Fokker-Planck) , it forms a closed system of partial differential equations which 
uniquely determines an associated Markovian diffusion process. The third 
Newton law in the mean is utilised to generate diffusion-type processes which 
are either anomalous (enhanced), or generically non-dispersive. 
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1 Smoluchowski diffusion processes, Stokes resis- 
tance and weakly out-of-equilibrium systems 

1.1 Traditional phenomenology 

In random media that at are statistically at rest, motion of single tracers or disper- 
sion of pollutants, also in the presence of external conservative force fields, may be 
consistently described in terms of Smoluchowski diffusion processes, |2|. Their for- 
ward drifts are interpreted to give account of the mean velocity j3~ 1 F/m imparted 
by the external conservative force to diffusing particles, an outcome of the accel- 
eration F/m "felt" and accumulated on the relaxation time scale This time 
scale is regarded to be much finer than the one appropriate for the coarse-grained 
description. The latter, At ^> is still significantly smaller than the available 
phenomenological (observational) resolution. 

Our basic intuitions are rooted in the theory of a Brownian motion suitable 
for the description of colloidal particles floating in a liquid. However an issue of 
particular "causes of diffusion" is not that relevant, and the proper arena for the 
Brownian motion can be not only a viscous fluid, or a dilute gas, but any inter- 
acting many-particle system and even less specific random medium with a suitable 
microscopic behaviour. The problem of "how to make a heat bath" needs to be 
addressed, |2|, ||. In the present paper, we take for granted the validity (even if 
diminished to the status of an approximate theory or the conceptual playground) of 
the standard Kramers and Smoluchowski diffusion process scenarios, as the reference 
mathematical models of random transport in the equilibrium bath. 

If we consider a fluid in thermal equilibrium as the noise carrier, a kinetic theory 
viewpoint amounts to visualizing the constituent molecules that collide not only with 
each other but also with the tagged (colloidal) particle, so enforcing its observed 
erratic motion. Clearly, random molecular collisions both initiate and maintain an 
incessant irregular motion of the tagged particle. Once the particle is in motion, 
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we need to account for an additional statistical effect of molecular impacts on the 
actually moving particle. It is phenomenologically encoded in the Stokes resistance, 
coming from averaging over a molecular "rain" along some portion of the trajectory, 
which is proportional to the velocity of the particle. In the phase-space description, 
the particle velocity V(t) = v is a random quantity. The damping effect is not 
random. It stands for a statistically accumulated, passive response of the medium. 
Indeed (see e.g. Section 2.1 in below), — -pt locally averages the genuine (ignored 
fine-detailed) random dynamics on the relaxation time scale (3~ l . 

For a bath in equilibrium, the internal relationship between the above random 
and systematic (due to friction) effects of a generic noise on the Brownian particle 
motion is the subject of fluctuation-dissipation theorems. This extends to situations 
when the external driving forces additionally modify the particle dynamics, but 
do not modify the statistics of noise (!). This feature is generic to Smoluchowski 
diffusions. (Conventionally, the fluctuation-dissipation theorems are interpreted, 
0, to set a general relationship between the response of diffusig particles in an 
equilibrium bath to an external force, and the internal fluctuation of the bath in the 
absence of this disturbance.) 

By means of Einstein's fluctuation-dissipation theorem we are given the diffusion 
coefficient D = It characterises an intensity of the spatial (Wiener) noise in 
terms of the friction parameter (3, with k the Boltzmann constant, T the equilib- 
rium temperature of the bath, and m mass of a diffusing particle, [[jj]. A formal 
exploitation of the Stokes formula (derivable on the basis of pure kinetic arguments 
from the Boltzmann theory, || ) transfers to the Brownian realm a concept of a fric- 
tional deceleration, originally suited for a macroscopic spherical particle of radius 
a and mass m moving in a fluid with viscosity coefficient rj. In the new context, 
it is the mean property of motion and not a particular single Brownian particle- 
in-motion attribute. Nonetheless, that allows to establish the value of the friction 
constant (3 = 6iir]a/m, the result amenable to positive experimental verifications in 
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the classic studies of the Brownian motion, [|T], 0. 



1.2 Meaning of stochastic models 

An observable Brownian motion of colloidal particles in a fluid, |l|, [7|], when inter- 
preted in terms of random processes, involves a number of mathematical subtleties, 
like e.g. an inherent nondifferentiability of sample paths in velocity space that reap- 
pears on the spatial arena of Smoluchowski processes as well. 

Clearly, the phase-space stochastic process is rather crude approximation of re- 
ality, if compared with the kinetic theory reasoning based on the explicit input of 
particle collisions to the Boltzmann equation and the related kinetic theory of hy- 
drodynamic flows, ||. The Boltzmann equation can be interpreted in terms of a 
jump Markov process simulating collision events, and the diffusion-type Karmers 
equation may arise only in a suitable scaling limit, ||. This limit can be justified 
in case of grazing collisions, or as a way to include the effect of long range forces 
by treating their influence in a statistical way: they should generically produce a 
sequence of small and almost random changes in the tracer particle velocity, 0. 

The Smoluchowski approximation takes us further away from the kinetic theory 
intuitions by projecting the phase-space theory of random motions into its configu- 
ration space image, ||, which is a spatial Markovian diffusion process. 

To quantify the above picture one usually departs from the phase-space (Langevin 
and Kramers) description of fluctuation phenomena. Let us consider an Ito equa- 
tion (in its symbolic differential version) for infinitesimal increments of the velocity 
random variable, exhibiting the systematic frictional resistance: 




(1) 



where W(t) stands for the normalised Wiener process. One can easily infer, |10|, 



the corresponding second Kolmogorov (Fokker-Planck) equation 



d t p(v , v, t) = (3DA ff + pVff ■ [vp(v , v, t)} 



(2) 
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for the transition probability density of the time homogeneous process in the velocity 
space alone, [II]. |3| . 
In view of: 

*■ ') = ( 2rt r(i- e -^ ) 3/2 ^W i-U > (3 > 

the time interval j3~ l effectively accounts for an approach of the transition density 
to the equilibrium Maxwell distribution. 

Let us consider an instantaneous velocity V(t) = v, that has been achieved in 
the course of the random evolution (1) beginning from a certain V(0) = vq. We 
can evaluate a conditional expectation value (local mean with respect to the law of 
random displacements (3)) over all randomly accessible velocities V(t + At) = v' at 
a time t + At, At > 0. It determines the forward drift of the process: 

b(v, t) = Um Am [J tfp(v, J, At) d 3 v' - v\ = -(3v (4) 

and thus, [0, 0, [T^|, provides us with an information about the mean tendency of 
the dynamics on small (but not too small if compared to f3~ x ) time scales. Evidently, 
that mean tendency in case of Eq. (1) is to decelerate the velocity v at the Stokes 
rate p—i ■ 

Sample paths of the Wiener process are nondifferentiable and play the role of 
(velocity space) idealizations of "true" trajectories. We need a conditional averaging 
over all paths emanating from v = V(t) to bypass, cf. Eq.(4), the nondifferentiability 
problem and then introduce the local deceleration rate (v — > v — /3vAt), ||. Clearly, 
the friction term and thus the deceleration concept, in the framework set by Eq. (1), 
arise as statistical, local mean value, quantities. All that must be sharply contrasted 
with the standard hydrodynamical meaning of the Stokes resistance formula, which 
refers to a genuine (the reduced Navier-Stokes equation is involved) fluid velocity 
field around a single uniformly propagating particle. It is the perturbation by the 
moving macroscopic body that gives rise to the force, with which the viscous fluid 
acts upon the concrete particle-like object. 



5 



Quite apart from the time-scale (At versus j3~ l ) issue, an individual Brownian 
particle neither moves uniformly, nor ever permanently relaxes to the state of rest. 
Actually, for a single Brownian particle, the respective velocity gain/loss (accelera- 
tion/deceleration rate) on the time scale is a random quantity. 

On the contrary, the standard interpretation of the Stokes resistance for a large 
particle in a viscous fluid, makes inavoidable an ultimate stopping of the particle, 
unless an external force would balance the damping and so maintain the (uniform) 
particle motion at a certain terminal velocity. This latter concept forms the basis 
of the Einstein fluctuation-dissipation formula in case of the gravitationally induced 
sedimentation phenomenon, cf. Jl|, albeit in case of the Brownian motion only local 
mean velocity fields (hence, ensemble averages) can be employed for that purpose. 



The picture of " a Brownian particle moving at its terminal velocity" , |I3] , is certainly 
inappropriate. 

— * — # 

We can supplement Eq. (1) by the spatial increment definition: dX(t) = V(t)dt, 
extending (1) to the phase-space process whose Smoluchowski projection in the large 
damping regime reduces to the pure spatial Wiener noise dX(t) = \/2DdW{t). In 
this case, the generic spatial scale (wandering distance over which the dissipation 
process in the mean is completed) is set by (D(3~ l ) l l 2 . The Smoluchowski ap- 
proximation amounts to ignoring the fine details of the dynamics on the and 
(-D/3 -1 ) 1 / 2 scales. One disregards all possible remnants of the phase space process 
that extend beyond those scales (the damping typically induces exp(—f3t) factors in 
all relevant formulas, hence a certain amount of "memory" must be eliminated to 
yield a Markov process), 

That particular disregard /ignore issue is worth emphasizing in the context of 
typical approximate (analytic or numerical) manipulations with the Brownian mo- 
tion. A strong solution of Eq. (1) (we take V^O) = v as the initial condition) has 
the form : 

V(t) =v exp(-(3t)+ [ exp[-(3(t- s)]V2DdW(s) . (5) 

Jo 
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This expression, albeit looking physically unrealistic, is a fairly accurate approx- 
imation of a phenomenological molecular collision scenario for the Brownian motion. 
Namely, let us introduce a finite Riemann sum approximation of the integral in Eq. 
(5): 

ft _ 

/ exp[-(t - s)]V2DdW(s) « exp(-/3t) V exp(J3nAt)AV n (6) 
Jo 

where AV n = \r2DAW(nAt) stands for the n-th consecutive velocity increment, 
e.g. an effect of all "random acceleration" events taking place in the time inter- 
val (nAt, (n + l)At). If we resort to a molecular collision mechanism (in a dilute 
gas for example), the velocity increments due to collisions of the Brownian particle 
with molecules of the bath can be viewed as statistically independent, and occur at 
an approximate rate 10~ 21 per second. The collisions may be interpreted to occur 
consecutively one after another, but multiple collision events are allowed as well. 
If we regard the time coarse-graining in Eq. (6) to refer to a typical relaxation 
time ~ i0~ 8 s, (notice that the respective coarse-graining appropriate for the 
Smoluchowski approximation, would involve At ^> it is obvious that each 

n-th velocity increment can be interpreted as a sum of an enormous number of inde- 
pendent identically distributed random variables (minute acceleration/deceleration, 
e.g. collision events). Let us stress that we exploit here a defining property of the 
Wiener process; the velocity random variable (5) of the Ornstein-Uhlenbeck process 



(l)-(3) is known not to have independent random increments, [10]. 

The accumulated fluctuation irregularities, on the (3~ l time scale, are the primary 
reason of the erratic behaviour of the Brownian motion in the velocity space. Thus, 
an integration with respect to the normalised Wiener process W(t) in Eq. (5), 
quite satisfactorily (in a suitable scaling limit) models a cumulative outcome of 
phenomenologically motivated impulses (kicks exerted by the noise carrier) on the 
Brownian particle. 

Let us emphasize at this point, 0, that it is not correct to think that conspicuous 
jiggles in the Brownian trajectory are due to single kicks. (In mathematical terms 
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the situation is even worse, since the Wiener process has an unbounded variation 
on aritrary time intervals.) The realistic Brownian motion is unbelievably gentle, 
specifically if we refer to a heavy particle immersed in a gas of much lighter ones. 
Each collision (kick) has an entirely negligible effect on the velocity of the Brow- 
nian particle. It is only fluctuations in the accumulation of an enormous number 
of very slight changes in the particle velocity which give the trajectory its irregu- 
lar appearence, both in velocity and position pictorial representations of motion. 
Conequently, equation (5) makes sense as a useful approximation of physical phe- 
nomena, if the coarse-grained time scale At in (6) (which is far below an obser- 
vational one) is much larger than /3~ l , but still small enough for a sufficiently fine 
coarse-graining (6) of (5). 

The previously mentioned space and time scales justify the utility of the Smolu- 
chowski description of conservative force effects upon a spatially diffusing particle, 

dX(t) = ^—dt + V2DdW(t) . (7) 

mp 

The Smoluchowski forward drift can be traced back to a presumed selective 
action of the force on the Brownian particle that has a negligible effect on the 
thermal bath. Indeed, we take for granted that there is no physically relevant mean 
(induced) flow in the bath proper, unless the isothermality assumption is abandoned, 
]T3| , [TJJ, or other hitherto disregarded effects on the (3~ l scale (like those due to the 
action-reaction principle, [|ll|]) are incorporated into the formalism. 

Brownian particles appear to be drifting on the local average relative to the bath, 
with a uniform (but in the mean !) velocity Clearly, a repeated series of observa- 
tions, at regularly spaced time intervals, of a single tagged Brownian particle would 
not in general reveal any specific motion tendency. The graphical picture of motion 
would be as irregular as ever ("no purpose" sample paths of the Wiener process, 
like in the absence of any force). The respective coarse-grained approximation of 
the trajectory does correspond to an experimental resolution, which is incomparably 
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rougher than the previous time scales. 

Only a numerical simulation of the statistical ensemble of sample paths with 
a controlled probability distribution (frequency, in fact) of initial conditions, or a 
realistic monitoring of the Brownian motion-induced dispersion of a low density 
cloud of dust particles suspended in a liquid, would lend a definite meaning to 
the "motion tendency" concept and to the related mean Brownian flows. Realistic 
diffusion processes (and diffusive transport) are observed in case of non-uniform 
concentrations of colloidal particles but they are regarded as a result of random 
migration of individual particles which is actually observable (under a microscope) 
and interpreted as the Brownian motion. 

In the context of a sedimentation phenomenon, a sufficiently long overall ob- 
servation time of a single tagged particle travelling along its erratic path (hours or 
days in a real, or computer simulation experiment), would presumably reveal that 
a particle more frequently visits certain spatial areas, in accordance with the baro- 
metric formula (ergodic features of motion). Early experiments on this issue, [TJ, [7|, 
pertained to a cloud of suspended particles executing an extremely slow (practi- 
cally adiabatic) diffusion process. In some cases, the spatial displacement of the 
size 0.2 ■ 10~ 6 m tracer particles has been measured in 30s time intervals, @, 0, to 
yield the observationally relevant outcomes. Very recent observations of an indi- 
vidual Brownian particle motion refer to 1 — 2.5 ■ 10 _6 m tracers and the time 
resolution. 

1.3 Problems with thermal equilibrium 

A spectacular solution of the sedimentation problem due to Smoluchowski, refers to 
the isothermal dynamics (1) which is constrained to the positive vertical semiaxis: 
dZ(t) = —(3cdt + P\ / 2DdW(t), with a reflecting boundary at the spatial point 0. 
The transition probability density of this one-dimensional process reads, JT], |T0|j : 

P{z,z ,t) = —-^—^{exp[-(z-z ) 2 /4Dt} + exp[-(z + z ) 2 /4Dt}}+ 
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C f°° 

-—=exp(-cz/D) J z+z ^ ct exp(-x 2 )dx (8) 

V 2(Dt) 1 / 2 

and shows that a particle starting its motion from any positive zo near 0, may 
wander along the positive semiaxis indefinitely. In particular, it can be transported 
against the gravitational force to an arbitrary height. Depending on the actual Zq 
the Brownian particle may have a higher probability to ascend than to descend. 
The corresponding mean work has been evaluated to be kT per particle, cf. [I], 10 



for a related discussion of the entropy decrease issue. Certainly, an ability of the 
medium to perform work (in the mean) i. e. to give a kinetic energy to the Brownian 
particle (on the ensemble average again), is not special to the sedimentation problem. 
It appears to be a universal feature of the thermal bath even in the absence of any 
external forces. 

In the equilibrium situation (t — > oo =>- p(z,z ,t) — > p(z)), we would arrive at 
the familiar balance condition: the mean tendency of motion (forward drift) due to 
gravitational acceleration must be exactly balanced by the oppositely directed mo- 
tion tendency of the diffusive (osmotic pressure, |l| [|) origin, valid for non- uniform 
concentrations of a contaminant in a solvent. The latter tendency involves sending 
particles away from the areas of higher probability of their presence (concentration, 
if a low density pollutant is considered) in accordance with the Fick formula for the 
diffusion current — c = —D vp Sv . The barometric formula does follow. 

It is interesting to spend a while on some tagged (single) particle features in 
the nonequilibrium-but isothermal-regime admitted by (8) for not too large times, 
when the Einstein (mean) balance condition is yet invalid. Clearly, to have defined 
an analogue of the Fick diffusion flow, a probability density of initial data zq must be 
chosen. In a computer simulation we would have under control a bunch of relevant 
sample paths (all consecutively executed on a fixed duration time interval [0,t]) 
and the related probability density evolution along the bunch. Depending on the 
initial data distribution, for a time period the osmotic drift would dominate the 
gravitationally induced drift. To this end, work (in the mean) must have been done 



10 



on Brownian particles at the expense of the bath. Local conversion of work into heat 
seems to be inavoidable, and local heat flows are inavoidable as well, cf. |jT3| , [7|. 
That creates a number of problems to which no attention is normally paid in the 
literature. 

Smoluchowski diffusions are conventionally regarded as isothermal processes (pos- 
sible heat flows are ignored by various reasons) . If we however admit that the emerg- 
ing tiny heat flows may have an affect on the particle transport due to the Brownian 
motion, a suitable description of the thermal inhomogeneities and their effects on 
the dispersion of Brownian particles must be invented. 

For example, in case of thermally inhomogeneous gases it is well known that a 
dust-free region appears about a hot body, showing that a temperature gradient 
has an effect on the motion of dust particles in a gas, ||16|| . Tracer particles are 
transported away from the hot areas to the cooler ones (we may interpret that as 
a repulsion by the heated domain). On the reverse, they appear to be attracted 
by the cooler areas while escaping from the hot domains. Particles float down the 
temperature gradients. 



Remark 1: Let us mention an approach to quite similiar problem, |L3[ which 
was originally formulated for a cloud of contaminants in a liquid, under the following 
phenomenological assumption: " a gas of Brownian particles falling in gravity should 
leave a trail of warm fluid in its wake, since its potential energy is being converted 
into heat" . Obviously, if the particles would move against gravitational force, then 
the temperature of the medium should locally drop down. Those features, if we 
are to keep track of the local heating and cooling (as opposed to the isothermal 
Einstein or Smoluchowski diffusive dynamics) were interpreted as a source of the 
space-time dependence of temperature. The Fokker-Planck equation must then be 
supplemented by an evolution equation for the temperature field, (a clear-cut kinetic 
theory reasoning can be read out in this strategy), so that the coupled nonlinear 
system would take the form of a "consistent thermodynamical system" i. e. the 
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one manifestly respecting the first and second laws of thermodynamics, see e. g. 
Ref. Here, the heat flows are assumed to be neither slow nor fast enough to be 
effectively disregarded. 



Remark 2: An issue of suitable slow and fast process time scales is crucial in 
our discussion. One should be aware that local temperature perturbations of the 
bath may be safely neglected when the dissipation (fast process) time scale and the 
diffusion (observational for tagged particles) time scales are generically incompatible, 
like e.g. in case of a rapid dissipation set against a slow diffusion process. Then, 
the usual isothermal diffusion process follows and the standard Brownian motion 
paradigm is left intact. 

In the discussion of the Carnot principle, in reference to the free Brownian motion 
and to the sedimentation problem, [l], [7j , Brownian particle fluctuations are regarded 
to occur due to causes that are intrinsic to the random medium. If we think about 
minute acceleration/deceleration events that modify (say, at a rate of 10 21 times 
per second) velocities of realistic particles, the microscopic energy-momentum con- 
servation laws need to be respected in each separate collision event. In contrast to 
derivations based on the Boltzmann colllision scenario, this feature is completely 
alien to the Brownian motion theory, cf. [II], O. (This happens quite apart from 
the elusive power of the flutuation-dissipation mechanism, [||: "The friction, or 
more generally the resistance of a given system, represents the method by which the 
external work is dissipated into microscopic thermal energy. The reverse process is 
the generation of random force as the result of thermal fluctuation".) 

Let us point out that the phenomenology allowing to regard Eq. (1) as a satis- 
factory model for Brownian particle velocity fluctuations, blurs possible advantages 
of the closely related argument appearing in the diffusion modeling of the Raileyigh 
piston, Hl7| : " Collisions between the large particle (piston) and the bath particles 
are elastic. After each collision the bath particles are given a new distribution, which 
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evolves until the next collision". 

To avoid an apparent contradiction with the law of energy conservation, there 
seems tempting to require that each minute acceleration of a Brownian particle is 
accompanied by a minute cooling of the medium in its immediate neighbourhood. 
Correspondingly, any deceleration event should induce a local heating of the im- 
mediate neighbourhood of a particle, see e.g. @. Since (cf. Eq. (1)), we always 
disregard the fine details of about 10 13 collision impacts on the Brownian particle 
on a typical relaxation time scale of 10 _8 s, there is definitely enough room to allow 
for local statistical measures of heating and cooling. 

Remark 3: The heating, cooling and temperature notions are quantitative mean 
measures of the degree of agitation of the noise carrier. If the random medium is 
interpreted on the molecular level to be composed of light particles, these measures 
can be related with the mean square deviation of the (bath molecule) velocity ran- 
dom variable, or an average of the squared velocity if its mean value vanishes. Those 
quanties are purely statistical characteristics of the bath and not of the Brownian 
particle immersed in it. Only under restrictive thermal equilibrium conditions, the 
notion of temperature appropriate to the bath, can be elevated to the status of a 
measure of a thermal agitation for tracer particles. 

1.4 Goals 

In the Kramers approach to the phase-space dynamics, the stochastic properties 
of the medium were considered to be independent from random phase-space data 
of the Brownian particle. The statistics of noise (e. g. the thermal equilibrium 
features of the bath) must have remained unperturbed by the very presence of the 
particle and its phase-space fluctuations, albeit those are enforced by the intrinsic 
randomness of the bath. One assumes that there is no relevant dynamical response 
of the bath to the very presence of the Brownian particle and its induced dynamics. 
We recall that the fluctuation-dissipation theorems merely account for a passive 
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response, in terms of the statistically implemented Stokes resistance of the bath to 
the particle-in-motion. 

On the other hand, the above local heating and cooling notions refer to the dy- 
namical response of the bath to the Brownian particle which is immersed in it. We 
need to account for the out-of-equilibrium properties of the bath in the presence 
of a single particle, whose motion is solely enforced by the bath. Thus very weak 
(although possibly fast) heat flows should accompany an individual Brownian par- 
ticle motion. That conforms with an obvious intuition that no physical system is 
ever in thermodynamical equilibrium, and normally it is necessary to idealize the 
situation by regarding "fast processes" to be completed, while the "slow" ones are 
still running. 

In such a weakly non-equilibrium system with small heat flows, we may expect 
that the standard equilibrium temperature notion is replaced by an effective tem- 
perature notion (and an effective thermal equilibrium), which depends on the chosen 
fast-versus-slow-process time scales and the ensemble averaging. Repetitions of a 
single particle experiment in the same thermal bath, should now be replaced by rep- 
etitions of the same experiment with different realizations of the out-of-equilibrium 
heat bath, cf. |18|] for similiar concepts in the context of randomly disordered media. 

In the above sense only an effective isothermal regime may be maintained, since 
the compensating heat flows allow to regain the equilibrium temperature almost 
instantaneously. After averaging over the tracer particle ensemble and the corre- 
sponding (weakly out-of-equilibrium) sample paths, we should be able to capture 
possible statistically relevant effects due to temperature inhomogoneities (small de- 
viations from thermal equilibrium conditions) and the resultant effective flows in a 
bath, that are possibly induced (via back reaction on the time scale probably larger 
than (3~ l but still much below the Smoluchowski approximation time scale At) by 
propagating Brownian particles. All that is to happen well beyond the limits of an 
available, cf. J?], observational resolution. 
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We have thus set a phenomenology of a specific dynamical response of the random 
medium to the Brownian motion of a single tagged particle, whose sole outcome are 
the minute deviations from the thermal equilibrium of the fluctuating medium itself 
and the resultant heat flows (needed to restore the equilibrium). This is independent 
from the traditional frictional resistance argument, directly referring to the thermal 
equilibrium conditions and the (3~ l time scale. 

We deal here with a generic feedback mechanism scenario. The Brownian particle 
propagates "at the expense" of the bath, which however remains "close" to its ther- 
mal equilibrium. The bath in turn actively reacts back to what is being happenning 
to the particle in the course of its propagation. The instantaneous local deviations 
from the state of equilibrium ("perturbation of noise") along the trajectory surely 
have an effect on each subsequent stage of the particle propagation. Even if residual 
for an individual Brownian particle and its sample path, the feedback effect is ex- 
pected to accumulate statistically (on the ensemble average) to a sizable quantity. 
Under the name of the "Brownian recoil principle" we have made a preliminary 
study of such random dynamics in Refs. O, |T2] . 

Remark 4: In connection with the previous Remark 1, let us mention that a 
spatial diffusion (Smoluchowski) approximation of the phase-space process, allows 
to reduce the number of independent local conservation laws (cf. ]15[ [L^, |20fl ) to 
two only. Therefore the Fokker-Planck (or continuity) equation can always be sup- 
plemented by another (independent) partial differential equation to form a closed 
system. Non-isothermal flow description needs to accomodate the variations of tem- 
perature of the bath, (cf. [p~3| , WW), while we investigate the limits of validity of 
the isothermal scenario. That amounts to inequivalent choices of the supplementary 
equation. We emphasize a single tagged particle-in a bath description in a repeat- 
able experimentation sequence, under basically the same (or very similiar) physical 
conditions. In such situation, a stochastic process executed by a single particle is 
practically isothermal ("the Brownian motion is unbelievably gentle", ||). Each 
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sample path of the Brownian particle (and the related cooling/heating phenomena 
induced along the path) is a random quantity. The state of the medium, if giving 
account of its dynamical response (i. e. deviation from thermal equilibrium) to the 
Brownian propagation, is a random quantity as well. It is the local mean (related to 
the ensemble average over various realisations of out-of-equilibrium conditions for 
the bath) that may properly quantify this picture. 

2 Local conservation laws for the Brownian mo- 
tion in the Smoluchowski approximation: diffu- 
sion currents and driving flows 

Before, we have identified the forward drift F — —VV as a quantitative measure 
of a statistical (local) tendency of the Brownian motion, obtained through averaging 
over an ensemble of sample paths. If we assign a probability density po(x) with which 
the initial data x = X(0) for Eq. (7) are distributed (weak solutions of (7) enter the 
scene), then the emergent Fick law would reveal a statistical tendency of particles 
to flow away from higher probability residence areas. This feature is encoded in the 
corresponding Fokker-Planck equation ( equivalent ly, a continuity equation): 

— * — » 

d t p = -V • (vp) = -V • [{— R - D^)p] (9) 
where a diffusion current velocity is v(x,t) = b(x,t) — D v ^_f '9 while the forward 



drift reads b(x, t) = cf. Eq. (7). Clearly, the local diffusion current (a local 



p(x,t) 

rrt/3' 

flow that might be experimentally observed for a cloud of suspended particles in 
a liquid) j = vp is nonzero in the nonequilibrium situation and a non-negligible 
matter transport occurs as a consequence of the Brownian motion, on the ensemble 
average. 

It is interesting to notice that the local velocity field v(x, t) obeys the natural 
(local) conservation law, which we quite intentionally pattern after the moment 
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identities (hierarchy of conservation laws) valid for the Boltzmann and Kramers 
equations, 0, [OJ. The pertinent momentum conservation law directly originates 
from the rules of the Ito calculus for Markovian diffusion processes, [fj|, and from 
the first moment equation in the diffusion approximation (!) of the Kramers theory, 
I HI: 

d t v + (v- V)v = V(fi-Q) . (10) 



drift b(x) = as follows: 



An effective potential function Q(x) can be expressed in terms of the forward 

n = £m + ^- p - (11 > 



Let us emphasize that it is the diffusion (Smoluchowski) approximation, |)], [15 
which makes the right-hand-side of Eq. (10) substantially different from the usual 
moment equations appropriate for the Brownian motion, ||15|| . In particular, the 
force F presumed to act upon an individual particle, does not give rise in Eq. (10) 
to the expression — ^VV which might be expected on the basis of kinetic theory 
intuitions and moment identities directly derivable from the Karmers equation, but 
to the term cf. Eq. (11). 

Moreover, instead of the standard pressure term, there appears a contribution 
from a probability density p-dependent potential Q(x,t). It is is given in terms of 
the so-called osmotic velocity field u(x,t), (cf. 0): 

Q(x,t) = ~u 2 + DV-u (12) 

u(x, t) = DV lnp(x, t) 

and is generic to a local momentum conservation law respected by isothermal Marko- 
vian diffusion processes, cf . 0, |11], |12|, [L9| . Notice that in case of the free Brownian 
motion (admitted, if we set Q = 0), we would have v(x,t) = —u(x,t) for all times. 

An equivalent form of the potential (12): Q = 2D 2 A ^/ 2 induces rather obvious 
quantum mechanical associations (the de Broglie-Bohm "quantum potential" with 
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the opposite sign, modulo an adjustment of constants), pif . In the context of the 
Brownian motion, this "quantum potential" has been deduced in earlier investiga- 
tions of local conservation laws, 



Remark 5: Let us notice that by demanding Q = Q identically for all x, t, 
we would reduce Eq. (10) to d t v + {y ■ V)v = 0. Despite of its classical-looking 
Riemann equation form, this conservation law still refers to a diffusion process. 
Namely, in view of (11), we must identify forward drifts with osmotic velocity fields 
and DWln p(x) = ^ holds true. The related diffusion process is stationary and 
preserves the probability measure (p is now time-independent). Cf. for more 
general considerations on that issue. 

As repeatedly stated before, Smoluchowski drifts refer to mean motions relative 
to the bath at rest, and there is no place for any flows intrinsic to the random 
medium in this formalism. On the other hand, it is of fundamental importance to 
understand how genuine flows in a random medium may be generated and what 



would be their effect on dispersion, [12\. (Solutions of the incompressible Navier- 
Stokes equation may serve as a common-sense model of the flow in a bath, and the 
diffusion enhancement is known to be related to various turbulent motion scenarios.) 

To analyze random perturbations that are either superimposed upon or are in- 
trinsic to a driving deterministic motion, a configuration space equation x = w(x, t) 
is normally invoked, which is next replaced by a formal infinitesimal representation 
of an Ito diffusion process 

dX(t) = b(X(t),t)dt + V2DdW(t) (13) 

patterned after Eq. (7). The tacit assumption (basically wrong, |p3[| ) is that b does 
not substantially differ from w. 

It is useful to exploit a standard phase-space argument that is valid, under 
isothermal conditions, for a Markovian diffusion process taking place in (or relative 
to) a flow w(x,t) with as yet unspecified dynamics nor concrete physical origin. We 
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account for an explicit force (here, acceleration K = F/m) exerted upon diffusing 
particles, while not directly affecting the driving flow itself. Namely, |TD|, |2J, let us 
set for infinitesimal increments of phase-space random variables: 

dX(t) = V(t)dt 

dV{t) = (3[w(x } t) - V(t)]dt + K(x)dt + (3V2DdW(t) . (14) 

Following the leading idea of the Smoluchowski approximation, we assume that 
P is large, and consider the process on time scales significantly exceeding (that 
is normally achieved by taking (3 to be very large, cf. the infinite friction limit 
procedure). Then, an appropriate choice of the velocity field w(x, t) may in principle 
guarantee, [Q], the convergence of the spatial part X(t) of the process to the Ito 
diffusion process with infinitesimal increments : 

dX(t) = [w(x, t) + -xK\dt + V2DdW{t) . (15) 

P 

Consequently, the forward drift of the process would read b(x,t) = w(x,t) + 
-sK{x). Notice that the fi~ x K contribution can be safely ignored if we are interested 
in the dominant driving motion. 

Throughout the paper we are interested in Markovian diffusion processes, which 
propagate respectively the phase-space or configuration space probability densities 
(weak solutions of stochastic differential equations are thus involved). In the con- 
figuration space variant corresponding to Eqs. (13), (15), we deal with a stochastic 
process whose probability density p(x, t) evolves according to the standard Fokker- 
Planck equation 

d t p = DAp - V • (bp) (16) 

which is supplemented by the momentum conservation law (in the mean) of the 
form (10) for v = b — D^£. If we compare Eq. (15) with Eq. (7), we realize 
that the transformation of drifts has been executed. Under suitable restrictions, we 
can relate probability measures corresponding to those equations by means of the 
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Cameron-Martin-Girsanov theory of measure transformations. The Radon-Nikodym 
derivative of measures is here involved and for suitable forward drifts that are gradi- 
ent fields it yields, [0, the most general form of an auxiliary potential Q(x, t) that 
is allowed to appear in Eq. (10): 

fi(f,t) = 2D[9 t 0+^ + V-6)] . (17) 

Here b(x,t) = 2DV<p(x,t). 

Eqs. (17) and (11) are trivial identities, if we take for granted that all drifts are 
known from the beginning, like in case of typical Smoluchowski diffusions where the 
external force F is a priori postulated. We can proceed otherwise and, on the con- 
trary, one can depart from a suitably chosen space-time dependent function fl(x,t). 
Then Eq. (16) should be considered as a nonlinear (Riccatti-type) equation which 
is to be solved with respect to the drift field potential $(a?, t). Such solution, when 
inserted to the Fokker-Planck equation (16) would ultimately yield an evolution of 
an initial probability density p(x, 0). 

From this point of view, while developing the formalism, one should decide what 
is a quantity of a primary physical interest: the field of drifts b(x, t) or the potential 
Q(x,t). They are not independent quantities, and enter the discussion as entangled 
objects. Mathematical features of the formalism appear to depend crucially on 
the properties (like continuity, local and global boundedness, Rellich class) of the 



potential Q, see e.g. [12 



If we decide that the momentum conservation law is governed by a bounded 



from below, continuous function Q(x,t), cf. [12]), then it seems worthwhile to 
mention a close connection of the considered framework with the general theory of 
small random perturbations of the classical Hamilton- Jacobi dynamics, [|19|]. An 
assumption that the forward drift is defined in terms of a gradient of a suitable 
function allows to rewrite the formula (17) in the form clearly reminiscent of the 
Hamilton- Jacobi equation (we set $ = 2D<p in Eq. (17) and take 0) as the 
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initial data for the t > evolution): 

ft = d t $ + i|V$| 2 + £>A$ . (18) 

An associated function (known as the so-called backward drift of a Markovian dif- 
fusion process) 6* = b — 2u, cf. Eq. (11), if denoted 6* = V$* is known to yield 



another modified Hamilton- Jacobi equation, |12| , |19 

n = d t $* + i|v$,| 2 -DA$, (19) 

to be solved with given terminal data &*(x, T) for times < t < T . 

Eq. (19) is identifiable as the so-called Hamilton- Jacobi-Bellmann programming 



equation in the optimal control of stochastic diffusion processes, [T{|. A related 
issue of viscosity solutions of the standard Hamilton- Jacobi equation has been ex- 
tensively studied in the literature as the D [ limit of solutions of the modified 
(e.g. Bellmann) equation. A direct connection (via the logarithmic Hopf-Cole trans- 
formation) of Eq. (19) with the forced Burgers equation and the generalised heat 
equation (hence, with the related Feynman-Kac potentials, semigroups and kernels) 



is well known, |12], [19], [23| . 

There is however more interesting to observe that a gradient field ansatz for the 
diffusion current velocity (v = VS): 

d t p = -V • l(VS)p] (20) 

allows to transform the momentum conservation law (10) of a Markovian diffusion 
process to the universal Hamilton- Jacobi form: 

n = d t s + ^\vs\ 2 + Q (21) 

where Q(x, t) was defined before in Eq. (12). By applying the gradient operation to 
Eq. (20) we recover (10). Notice that Eq. (20) is sensitive to any additive (constant 
or time- dependent) modification of the potential Q. In the above, the contribution 
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due to Q is a direct consequence of an initial probability measure choice for the 
diffusion process, while Q via Eq. (17) alone does account for an appropriate forward 
drift of the process. 

The derivation of a hierarchy of local conservation laws (moment equations) for 
the Kramers equation can be patterned after the standard procedure for the Boltz- 
mann equation, || || [L5[]. Those laws do not form a closed system and additional 
specifications (like the familiar thermodynamical equation of state) are needed to 
that end. In case of the isothermal Brownian motion, when considered in the large 
friction regime (e.g. Smoluchowski diffusion approximation), the Fokker-Planck 
equation must be supplemented by one conservation law only to yield a closed sys- 
tem. Such system uniquely determines the stochastic process. 

That happens under a definite choice of external forces, and hence Smoluchowski 
drifts. If the drifts are not a priori specified, then the only freedom left in the 
momentum conservation law amounts to the choice of a concrete functional form 
for the potential Q(x,t). In the theory of Brownian motion this particular decision 
making replaces the standard equation of state constraint, suitable for the kinetic 
theory description of gases and liquids. 

In view of more sensitive dependence on the potential and hence more detailed 
discrimination between distinct dynamics scenarios, we adopt the Hamilton- Jacobi 
equation (21) as a generic substitute of the momentum conservation law (10). Thus, 
we can consider a closed system which is composed of the continuity equation d t p = 
— W(vp) (this, in view of v = b — D^, is equivalent to the Fokker-Planck equation 
(16)) and the Hamilton- Jacobi equation (21), plus suitable initial (and/or boundary) 
data. Alternatively, we can supplement the Fokker-Planck equation (16) by the 
nonlinear Riccatti-type equation (17) to form a closed system of partial differential 
equations, provided the functional form of Q is pre-selected. In contrast to the pair 
(16), (21) where p(x,t) enters an entangled relationship, the pair (16), (17) is not 
entangled. 
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We need to stress that it is the closed system of Eqs. (20) and (21) which directly 
refers to physically motivated local conservations laws (moment equations) associ- 
ated with the Brownian motion, PO, [TR O , and to the respective diffusion currents. 
The underlying Markovian diffusion process is then specified uniquely (that would 
not be the case if (10) has been used instead of (21)). However, Eqs. (20) and (21) 
form a coupled nonlinear system, whose analytic solutions are not readily accessi- 
ble. It is therefore useful to know that a linearisation of this formidable nonlinear 
problem is provided by a time-adjoint pair of generalised diffusion equations (cf. 
the Appendix) in the framework of the Schrodinger boundary data problem. The 
standard heat equation appears as a very special case in this formalism. 



3 The third Newton law in the mean 

3.1 Free Brownian motion in terms of local conservation 
laws 

The local conservation law (10) acquires a direct physical meaning (the rate of change 
of momentum carried by a locally co-moving with the flow volume, only if av- 
eraged with respect to p(x, t) over a simply connected spatial area. Namely, if V 
stands for a volume enclosed by a two-dimensional outward oriented surface dV, 
we define a co-moving volume on small time scales, by deforming the boundary 
surface in accordance with the local current velocity field values. Namely, we con- 
sider at time t the displacement of the boundary surface dV(t) defined as follows: 
x G dV — > x + v(x,t)At for all x G dV. Up to the first order in At this guaran- 
tees the conservation of mass (probability measure) contained in V at time t i. e. 
Sv(t+At) P& t + At)d 3 x - J v{t) p(x, t)d 3 x ~ 0. 

The corresponding (to the leading order in At) quantitative momentum rate- 
of-change measure reads, cf. |T2[], J v pV(Q — Q)d 3 x. In view of diQ = ^J2jdjPij, 
where the stress tensor = D 2 pdidj In p is determined up to an additive time- 
dependent or constant term, the standard divergence theorem allows to isolate an 
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explicit surface (due to stresses or of the pressure-type) contribution. Namely, there 
holds — J v pdiQd 3 x = — j dv J2j Pijdaj, with da being an infinitesimal area element 
of dV in R 3 . 

For a particular case of the free Brownian expansion: 

= J^W eXp{ -^ = [4ttD(*U)] 3 / 2 eXPl ~mTTo) ] (22) 

where a 2 = 4_Dt and t > 0, t > 0, we would have 

P y (f , f) = P(x , t)5 l3 = - ^JL_ p (x, t)5 i3 (23) 

for all x E R 3 and t > 0. Here % stands for the Kronecker symbol. Then, 
— J y pVQd 3 x = — J dv P da , where 

r 2 3D 
Q(f , t) = - - - , . (24) 

1 ^ 8(t + t ) 2 2(* + * ) 1 ; 

The current velocity v(x,t) = VS(x,t) = 2 (t+t ) ' a P ar ^ f rom solving <9 t p = 
— V ■ (-up) and <9 t -u + (#• V)v = —VQ with p and v = —u standing for initial data, 
is also linked to the Hamilton- Jacobi equation 

d t S + ^\VS\ 2 + Q = (25) 

whose solution is: S(x,t) = 4( - f + fo ) + ^Dln[AnD(t + to)]. 

Let us observe that the initial data vo = —DVln po = — uo for the current 
velocity field indicate that we have totally ignored a crucial preliminary stage of the 
dynamics on the /3" 1 time scale, when the Brownian expansion of an initially static 
ensemble has been ignited and so particles have been ultimately set in motion. 

Notice also that our "osmotic expansion pressure" P{x, t) is not positive definite, 
in contrast to the familiar kinetic theory (equation of state) expression for the pres- 
sure P(x) = a/(x), a > appropriate for gases. The admissibility of the negative 
sign of the "pressure" function encodes the fact that the Brownian evolving concen- 
tration of particles generically decompresses (blows up), instead of being compressed 
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by the surrounding medium. The compression (pressure upon the control volume 
coming from its surrounding) is the standard feature in the kinetic theory of gases, 
except for the cavitation phenomenon in liquids and the exotic blow-up conditions 
in the concentrations of dense hot matter. 

The loss (in view of the "osmotic" migration) of momentum stored in a con- 
trol volume at a given time, may be here interpreted in terms of an acceleration 
— J v pVQd 3 x induced by a fictituous "attractive force". By invoking an explicit 
Hamilton- Jacobi connection (21), we may attribute to a diffusing Brownian ensem- 
ble floating through a (locally co-moving) control volume V, the mean kinetic energy 
per unit of mass J v p^iPcPx. We can also evaluate the mean total kinetic energy 
per unit of mass obtained after extending integrations from V to R 3 . 

For the considered example, in view of < x 2 >= 6D(t+to), we have J R3 p\v 2 d 3 x = 
4 (t+ to y Notice that the mean energy f v pQ^v 2 + Q)d 3 x needs not to be positive. In- 
deed, this expression identically vanishes after extending integrations from V to R 3 . 
On the other hand the kinetic contribution, initially equal J R3 \pv 2 d 3 x = 3D /a 2 and 
evidently coming from nowhere, continually diminishes and is bound to disappear in 
the asymptotic t — > oo limit, when Brownian particles become uniformly distributed 
in space. 

3.2 The third Newton law in the mean and the Brownian 
recoil principle 

Normally, diffusion processes yielding a nontrivial matter transport (diffusion cur- 
rents) are observed for a non-uniform concentration of colloidal particles. We can 
devise a thought (numerical) experiment that gives rise to a corresponding trans- 
port in terms of an ensemble of sample Brownian motion realisations on a fixed 
finite time interval, instead of considering a multitude of them (migrating swarm 
of Brownian particles) simultaneously. One may surely implant particles at initial 
(random) space locations to mimic a certain probability density and next release 
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(individually in each sample experiment) and allow them to execute their Brownian 
paths independently, in a fixed duration time interval. In terms of such a particle 
ensemble, we can safely return back to the previous colloidal diffusion picture where 
migrating Brownian particles are also regarded as independent (and so are their 
individual Brownian motions). 

Consequently, in both visualizations, after the relaxation time the diffusion 
current is initiated " at the expense" of the bath. A nonzero mean kinetic energy must 
have been initially transferred (pumped) from the bath to the diffusing (blowing-up, 
expanding) swarm of particles. Only asymptotically this energy is being returned 
back to the bath. 

Recalling our previous discussion, for a sufficiently fast diffusion process, all that 
should correspond to a local cooling of the bath and implement a tiny deviation from 
its thermal equilibrium conditions in each single particle propagation (simulation) 
experiment. Accordingly, the tendency to regain the local thermal equilibrium by the 
bath (reflecting an attraction of tagged particles by the cooler areas) must result in 
induced local flows - they can become identifiable only on the ensemble average. All 
that is to happen not literally on the [3" 1 relaxation time scale but on (still relatively 
small) time scale At of the diffusion process which is well below the observational 
one. 

To each executed sample path there corresponds a sample realisation of the 
random medium (pushed slightly away from its thermal equilibrium, in view of 
the postulated feedback mechanism). Those random (sample) realisations of the 
bath should be ensemble-averaged as well to yield an "effective" bath in thermal 
equlibrium, which is however no longer in a statistical state of rest. The emergent 
driving flows mimic, on the average, the "return to equilibrium" of the bath in 
each sample propagation experiment. The thermal conditions are maintained on 
the ensemble average, so that the effective process guiding the ensemble dynamics 
can be viewed as isothermal. 
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We recall close links of this scenario with an idea of a random walk in a random 
medium, ||18|| . However, presently a disorder in the random medium is coupled to 
the randomness of the walk and is no longer of independent origin. 

Let us assume that "an effort" (hence, an energy loss) of the random medium, 
on the (3~ l scale, to produce a local Brownian diffusion current v(x, to) out of the 
initially static ensemble and thus to decompress (lower the blow-up tendency) an 
initial non-uniform probability distribution, results in the effective osmotic reaction 
of the random medium. Whatever is being transported away (on the ensemble 
average) according to the Fick law, is assumed to induce a compensating osmotic 
counterflow in the effective thermal bath. That is the Brownian recoil effect of Ref. 

n . 

Now, the particle swarm propagation scenario becomes entirely different from 
the standard one, (10), (20). First of all, the nonvanishing forward drift b = u 
is generated as a dynamical (effective, statistical here !) response of the bath to 
the enforced by the bath particle transport with the local velocity v = —u. Sec- 
ond, we need to account for a parellel inversion of the pressure effects (compression 
+VQ should replace the decompression — VQ) in the respective local momentum 
conservation law. 

Those features can be secured through an explicit realization of the action- 
reaction principle ("the Brownian recoil" effect, cf. Ref. |ll]]), which we promote to 
the status of the third Newton law in the mean. 

On the level of Eq. (10), once averaged over a finite volume, we interpret the 
momentum per unit of mass rate-of-change J v pV(Q — Q)d 3 x which occurs exclu- 
sively due to the Brownian expansion, to generate a counterbalancing rate-of-change 
tendency in the random medium. To account for the emerging forward drift and an 
obvious modification of the subsequent dynamics of an ensemble of (tagged) par- 
ticles, we re-define Eq. (10) by setting — J v pV(Q — Q)d 3 x in its right-hand-side 
instead of + J v pV(Q — Q)d 3 x . That amounts to an instantaneous realisation of 
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the third Newton law in the mean (action-reaction principle). Hence, the momen- 
tum conservation law for the process with a recoil (the reaction term replaces the 
decompressive "action" term) would read: 

d t v+(v- V)v = V(Q-Q) (26) 

implying that 

d t S + ^\VS\ 2 -Q = -n (27) 



stands for the corresponding Hamilton- Jacobi equation, cf. II9L E3, instead of Eq. 



(21). A suitable adjustment (re-setting) of the initial data is here necessary, which 
we shall explain in below (cf. Section 3.3). 

The major idea of the Brownian motion with a recoil is thus given. 

In the coarse-grained picture of motion we shall deal with a sequence of re- 
peatable scenarios realised on the Smoluchowski process time scale At: the Brown- 
ian swarm expansion build-up is accompanied by the parallel counterflow build-up, 
which in turn modifies the subsequent stage of the Brownian swarm migration (be- 
ing interpreted to modify the forward drift of the process) and the corresponding 
built-up anew counterflow. 

Although the new closed system of partial differential equations (20) and (27) is 
very different from the previous one (20), (21), nonetheless it describes Markovian 
diffusion-type processes again, |2|, |1^, ^5|. The link is particularly obvious if we 
observe that the new Hamilton- Jacobi equation (27) can be formally rewritten in 
the previous form (21) by introducing: 

tt r = d t S + hvS\ 2 + Q 

n r = 2Q-n (28) 

where Q represents the previously defined potential function of any Smoluchowski 
(or more general) diffusion process, (11). It is Q r which via Eq. (17) would determine 
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forward drifts of the Markovian diffusion process with a recoil. They must obey the 
Cameron-Martin-Girsanov identity 

n r = 2Q-n = 2D[d t <t> + i(J^ + V • 6)] . (29) 

Our system of equations (20), (27) is badly nonlinear and coupled, but its lineari- 
sation can be immediately given in terms of an adjoint pair of Schrodinger equations 
with a potential fl, [0, 0. Indeed, 

id t iP = —DAip + ^ (30) 



with a solution 



if) = p 1/2 exp{tS) (31) 



and its complex adjoint makes the job, if we regard p together with S to remain in 
conformity with the previous notations of Eqs. (20), (27). The choice of ip(x, 0) gives 
rise to a solvable Cauchy problem. This feature we shall exploit in below. Notice 
that, in view of Eq. (30), for time-indepedent Q, the total energy J r3 (^y — Q + pd 3 x 
of the diffusing ensemble is a conserved quantity. 

The problem (20), (27), (28) can be reformulated as the Schrodinger boundary 
data problem (cf. the Appendix), but the resulting generalised diffusion equations 
are nonlinearly coupled (by means of the potential Q r replacing the Feynman-Kac 
potential Q. Hence the previous linearisation needs to be exploited anyway. The 
general existence criterions for Markovian diffusion processes of that kind, were 
formulated in Ref. P5|, see also pi O and the Appendix. 



3.3 Brownian motion with a recoil as an anomalous (en- 
hanced) diffusion model 

For the clarity of discussion, we shall confine our considerations to one-dimensional 
problems. In the absence of external forces, we may consider a solution of equa- 
tions (in space dimension one) d t p = —V(vp) and d t v + (v V)t> = +VQ, where an 
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and 



initial probability density po(x) is chosen in correspondence with the previous free 
Brownian motion example. We denote a 2 = ADto- Then, 

OL x 2 c? 

P(X ' l) = i^+lOTp eXp[ ~a* + DH2 ] (32) 

,/ n f \ t x 2D(a 2 -2Dt)x 
b(x, t) = v(x, t) + u(x, t) = A _ (33) 

are the pertinent solutions. Notice that u(x, 0) = — = b(x, 0) amounts to 
f(x,0) = 0, while in the previous free Brownian case the initial current velocity 
was equal to —DVln p$. This re-adjustment of the initial data can be interpreted in 
terms of the counterbalancing (recoil) phenomenon: the would-be initial Brownian 
ensemble current velocity vq = ~uq is here completely saturated by the emerging 
forward drift bo = «o, see e.g. also [[□]]. This implies u(x,t) = — a l^^ t i and 
v(x,t) = a i^^l t i . Notice that VQ = j^^im , to be compared with the respec- 
tive expression — VQ = — a ?+£ Dt in the previous section. Presently, we deal with a 
fictituous "repulsive" force, which corresponds to the compression (pressure upon) 
of the Brownian ensemble due to the counter-reaction of the surrounding medium. 
We can write things more explicitly. Namely, now: 

and the corresponding pressure term (VQ = ^V P) reads 

giving a positive contribution +VQ to the local conservation law (26). 
The related Hamilton- Jacobi equation 

d t S+^\VS\ 2 = +Q (36) 

is solved by 

, 2D 2 xH „ . 2Dt. 

t] = a 4 + ADH 2 - arCtaU ( ~^ } (37) 
30 



With the above form of Q(x,t) one can readily check that equations (28) are 
identically satisfied, and that the Cameron-Martin-Girsanov constraint euqation for 
the forward drift of the Markovian diffusion process with a recoil is automatically 
valid for (ft = |Zn p + S: 

2Q = 2D[d t <j> + ±(-^ + V-b)} (38) 

cf. the general identity (29). 

In anology with our free Brownian motion discussion, let us observe that presently 

< x 2 >= — + . 39 

2 or 

It is easy to demonstrate (use a linearisation (30), (31) of the problem) that 
the quadratic dependence on time persists for arbitrarily shaped initial choices of 
the probability distribution po(x) > 0. That signalizes an anomalous behaviour 
(enhanced diffusion) of the pertinent Markovian process when Q = i. e. Q r = 2Q. 
We can evaluate the kinetic energy contribution 

ADH 2 



j p V —dx 



a 2 {a 4 + 4D 2 t 2 ) ^ 
which in contrast to the Brownian case shows up a continual growth up to the 
terminal (asymptotic) value This value was in turn an initial kinetic contribution 
in the previous Brownian example. In contrast to that case, the total energy integral 
is now finite (finite energy diffusions of Ref . |25| ) and reads 



2 



x \v 2 - Q)pdx = ^ (41) 
r 2 cr 

for all times (it is a conservation law). The asymptotic value of the current velocity 
v ~ | is twice larger than this appropriate for the Brownian motion, v ~ ^. 

3.4 Response to an external force: non-dispersive diffusion- 
type processes 

Let us regard F(x) = — W(x) as an external force field, whose effects on the 
dynamics of Brownian particles is encoded in the Smoluchowski diffusion process 
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equations (7), (9) and next (9)-(ll). Those in turn can be motivated by invoking 
the Kramers equation (14) and its Smoluchowski diffusion approximation (15). 

As emphasized before, on the level of local conservation laws, in the diffusion ap- 
proximation, the microscopic force F is represented by the (Feynman-Kac) potential 
Q, defined through the Girsanov formula (11). 

Let us adopt the third Newton law in the mean and the related Brownian recoil 
strategy (26)- (29) to this case. Evidently, the potential Q will explicitly appear in the 
linearization (30) of the problem. On the other hand, it is the potential Q r = 2Q — Q 
which via Eq. (29) determines forward drifts appropriate for the diffusion process 
with a recoil. In view of an inherent nonlinearity of the problem, one should not 
expect that the emergent drifts would allow for a simple decomposition met in Eqs. 
(14), (15). 

Our further discussion will be carried out in one space dimension and will focus 
on quadratic potentials. 

For a parabolic (harmonic oscillator) potential V(x) = ^moj 2 x 2 defining the 
acceleration K(x) = —uj 2 x, the corresponding Feynman-Kac potential (11) reads 
£l(x) = \ r y 2 x 2 — D'y, 7 = It is useful to mention that the choice of the repulsive 
potential V(x) = +\muj 2 x 2 would yield an innocent-looking modification by a con- 
stant in the function (17): Q(x) = \ r ) 2 x 2 + D^(. That demonstrates an extraodrinary 
sensitivity of the Riccatti-type equations (11) and (17) on the choice of Q. 

In fact, a suitable additive modification of |7 2 x 2 by a constant, allows to gener- 
ate (by directly solving the Riccatti-type equation (11)) the whole family of forward 
drifts pertaining to inequivalent stationary diffusion processes, cf. [0. Nonetheless, 
all of them correspond to the same VQ = +7 2 x generic contribution to the local 
momentum conservation law (10). Clearly, the law (10) does not provide a suffi- 
ciently fine discrimination between admissible stochastic motion scenarios, unless 
we know the Smoluchowski force and its potential from the beginning. It is only 
the Hamilton- Jacobi equation level, where the closed system of partial differential 
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equations (respectively (20), (21) or (20), (27)) determines the process uniquely. 

It is clear that stationary processes are the same both in case of the standard 
Brownian motion and the Brownian motion with a recoil. The respective propaga- 
tion scenarios substantially differ in the non-stationary case only. 

To exemplify the above statement (we have discussed before the Q = case), let 
us consider an explicit solution of Eqs. (20) and (28) in case of Q(x) = \^j 2 x 2 — -D7. 
By means of the linearisation (30), (31), this can be easily accomplished, cf. [p4| . 
We shall utilse exactly the same inital probability density po(x) as before. We have: 



p(x,t) 



W 11/2. 



7T 



'sin 2 (^2Dt) + 1 a i cos 2 (^2Dt)y 



' 6XP ^ sin 2 (^f2Dt) + 1 a' i cos 2 (^f2Dt^ ^ 

and 

s{xt) ^ DbflWxHg^Dt) , D^x 2 t 



(tan 2 (^2Dt + -fa 4 )sin(^2Dt) tan(^y2Dt) 
tan{^/j2Dt) 

The forward drift of the corresponding diffusion-type process reads: 



+Darctan [ — -] - D^t . (42) 



v ; sin 2 (^2Dt) +^a 4 cos 2 {^2Dt) v ; 

and X(t) = x holds true in terms of the random variable of the process. The 
additive decomposition of the drift, valid in Eqs (9), (15), is completely destroyed 
by the Brownian recoil scenario. Notice that b(x, 0) = — 2 -^f = u(x, 0), while (cf. (9)) 
b = -^5 = —"sx would hold true for all times, in case of the standard Smoluchowski 

rap 1 ' 

diffusion process. 

Because of the harmonic attraction and suitable initial probability measure 
choice, we have here wiped out all previously discussed enhanced diffusion features. 
Now, the dispersion is attentuated and actually the non-dispersive diffusion-type 
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process is realised: < x 2 > does not spread at all despite of the intrinsically stochas- 
tic nature of the dynamics (finite-energy diffusions of Ref . ) . 
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Appendix: Reconstruction of a Markovian diffusion process from the 
input-output statistics data 

There are many procedures to reproduce the intrinsic dynamics of a physical 
system from observable data, like e.g. the time series analysis. We shall outline an 
algorithm allowing to reconstruct the most likely microscopic motion scenario under 
an additional assumption that the sought for dynamics actually is a Markovian 
diffusion process. This reconstruction method is based on solving the so-called 
Schrodinger boundary-data and interpolation problem, [|19| , |T^ , [23f] . 

Given two strictly positive (usually on an open space-interval) boundary prob- 
ability densities Po(x), Pt(x) for a process with the time of duration T > 0. One 
can single out a unique Markovian diffusion process which is specified by solving the 
Schrodinger boundary data problem: 

m T (A,B)=f*xf*vm I {2,t) 

J A JB 

J d 3 xm T (x,y) = p T (y) 
where the joint probability distribution has a density 

m T (x, y) = u (x) k(x, 0, y, T) v T (y)0 
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and the two unknown functions Uo(x), Vt(v) come out as (unique) solutions, of the 
same sign, of the integral identities. To this end, we need to have at our disposal 
a continuous bounded strictly positive (ways to relax this assumption are known) 
function k(x, s, y, t), < s < t < T, which for our purposes (an obvious way to 
secure the Markov property) is chosen to be represented by familiar Feynman-Kac 
integral kernels of contractive dynamical semigroup operators: 

k(y,s,x,t) = J exp[- c(u(T),T)dT]dn^(u) 

The pertinent (interpolating) Markovian process can be ultimately determined 
by means of positive solutions (it is desirable to have them bounded) of the adjoint 
pair of generalised heat equations: 

d t u(x, t) = uAu(x, t) — c(x, t)u(x, t) 

d t v(x, t) = —uAv(x, t) + c(x, t)v(x, t) . 

Here, a function c(x, t) is restricted only by the positivity and continuity demand for 
the kernel. In the above, dpL% s l{u) is the conditional Wiener measure over sample 
paths of the standard Brownian motion. 

Solutions, upon suitable normalisation give rise to the Markovian diffusion pro- 
cess with the factorised probability density p(x, t) = u(x, t)v(x, i) which, while evolv- 
ing in time, interpolates between the boundary density data p(x, 0) and p(x, T). The 
interpolation admits an Ito realisation with the respective forward and backward 
drifts defined as follows: 

Vu(f, t) 



b(x,t) = 2u- 



v(x, t) 
Vu(x, t) 



K(x,t) = -2u , 

u(x, t) 

in the prescribed time interval [0, T]. 

For the forward interpolation, the familiar Fokker-Planck (second Kolmogorov) 
equation holds true: 

d t p(x,t) = uAp(x,t) — V[b(x,t)p(x,t)] 
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with p(x, 0) given, while for the backward interpolation (starting from p(x,T)) we 
have: 

d t p(x, t) = -uAp(x, t) - V[6*(x, t)p(x, t)] . 

— * 

The drifts are gradient fields, curlb = 0. As a consequence, those that are 
allowed by any prescribed choice of the function c(x, t) must fulfill the compatibility 
condition 

1 b 2 

c(x,t) = d t <5> + -(— + V6) 

which establishes the Girsanov-type connection of the forward drift b(x, t) = 2z/V<3>(x, t) 
with the Feynman-Kac potential c(x,t). In the considered Schrodinger's interpola- 
tion framework, the forward and backward drift fields are connected by the identity 

— * — * 

b* = b — 2vVlnp. 

For Markovian diffusion processes the notion of the backward transition proba- 
bility density p*(y, s, x, t) can be consistently introduced on each finite time interval, 
say < s < t < T: 

p(x, t)p*(y, s, x, t) = p(y, s, x, t)p(y, s) 

so that fp(y,s)p(y,s,x, t)d 3 y = p(x,t) and p(y,s) = J p*(y, s, x,t)p(x,t)d 3 x. 

The transport (density evolution) equations refer to processes running in opposite 
directions in a fixed, common for both time-duration period. The forward one 
executes an interpolation from the Borel set A to B, while the backward one executes 
an interpolation from B to A. 

The knowledge of the Feynman-Kac kernel implies that the transition probability 
density of the forward process reads: 

P{y, s , x, t) = k{y, s, x, t) V ^ X J - . 

v(y, s) 

while the corresponding transition probability density of the backward process has 
the form: 

( ^ _ v u ^ .My,s) 
P*{y,s,x,t) = k{y,s,x,t) 



u(x, t) 
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Obviously in the time interval < s < t < T there holds: 

u(x,t) = J u (y)k(y,s,x,t)d 3 y 

Consequently, we have fully determined the underlying (Markovian) random mo- 
tions, forward and backward, respectively. All that accounts for perturbations of 
(and conditioning upon) the Wiener noise. 



Remark 6: Various partial differential equations associated with Markovian 
diffusion processes are known not to be invariant under time reversal (hence being 
dissipative and linked to irreversible physical phenomena). However, the correspod- 
ing processes admit a statistical inversion. Let us consider a process running in a 
finite time interval, say [0, T\. We may consistently define a process running back- 
ward in time in this interval and reproducing the most likely (statistical) past of 
the process, given the present probability measure data. See e.g. |fj, [11], [12], [0J and 
26[ . In fact, cf. |27| p. 255: " any probabilistic treatement of the heat equation 



involves a time-reversal". This feature is explicitly utilized in the analysis of the 
above outlined Schrodinger boundary-data and interpolation problem, [19|, [L2|] . 
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